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Abstract 

In this paper, we extend the jump-diffusion model proposed by Davis and Lleo to include 
jumps in asset prices as well as valuation factors. The criterion, following earlier work by 
Bielecki, Pliska, Nagai and others, is risk-sensitive optimization (equivalent to maximizing 
the expected growth rate subject to a constraint on variance.) In this setting, the Hamilton- 
Jacobi-Bellman equation is a partial integro-differential PDE. The main result of the paper 
is to show that the value function of the control problem is the unique viscosity solution of 
the Hamilton- Jacobi-Bellman equation. 

Keywords: Asset management, risk-sensitive stochastic control, jump diffusion processes, 
Poisson point processes. Levy processes, HJB PDE, pohcy improvement. 

1 Introduction 

In this paper, we extend the jump diffusion risk-sensitive asset management model proposed by 
Davis and Lleo [19J to ahow jumps in both asset prices and factor levels. 

Risk-sensitive control generalizes classical stochastic control by parametrizing explicitly the 
degree of risk aversion or risk tolerance of the optimizing agent. In risk-sensitive control, the 
decision maker's objective is to select a control policy h{t) to maximize the criterion 

J{t,x,h;e) := -^InE [e-''^^*'^''^)] (1) 

where t is the time, x is the state variable, F is a given reward function, and the risk sensitivity 
9 g] — 1,0[U]0, cx)) is an exogenous parameter representing the decision maker's degree of risk 
aversion. A Taylor expansion of this criterion around 6 = yields 

J(t, X, h;9)=E [F{t, X, h)] - ^Var [F{t, x, h)] + 0{e^) (2) 

which shows that the risk-sensitive criterion amounts to maximizing E [F{t, x, h)] subject to a 
penalty for variance. Jacobson [28J, Whittle Bensoussan and Van Schuppen [S] led the 
theoretical development of risk sensitive control while Lefebvre and Montulet [3JJ , Fleming |25] 
and Bielecki and Pliska jllj pioneered the financial application of risk-sensitive control. In par- 
ticular, Bielecki and Pliska proposed the logarithm of the investor's wealth as a reward function, 
so that the investor's objective is to maximize the risk-sensitive (log) return of his/her portfolio 
or alternatively to maximize a function of the power utility (HARA) of terminal wealth. Bielecki 
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and Pliska brought an enormous contribution to the field by studying the economic properties of 
the risk-sensitive asset management criterion (see [l^), extending the asset management model 
into an intertemporal CAPM (^Mj), working on transaction costs ([Ijj), numerical methods 
(|10|) and considering factors driven by a CIR model (|15]). Other main contributors include 
Kuroda and Nagai [31] who introduced an elegant solution method based on a change of measure 
argument. Davis and Lleo applied this change of measure technique to solve a benchmarked 
investment problem in which an investor selects an asset allocation to outperform a given fi- 
nancial benchmark (see [18] ) and analyzed the link between optimal portfolios and fractional 
Kelly strategies (see [2D]). More recently, Davis and Lleo [12] extended the risk-sensitive asset 
management model by allowing jumps in asset prices. 

In this chapter, our contribution is to allow not only jumps in asset prices but also in the 
level of the underlying valuation factors. Once we introduce jumps in the factors, the Bellman 
equation becomes a nonlinear Partial Integro-Differential equation and an analytical or classical 
C^''^ solutions may not exist. As a result, to give a sense to the relation between the value func- 
tion and the risk sensitive Hamilton-Jacobi-Bellman Partial Integro Differential Equation (RS 
HJB PIDE), we consider a class of weak solutions called viscosity solutions, which have gained 
a widespread acceptance in control theory in recent years. The main results are a comparison 
theorem and the proof that the value function of the control problem under consideration is the 
unique continuous viscosity solution of the associated RS HJB PIDE. In particular, the proof 
of the comparison results uses non-standard arguments to circumvent difficulties linked to the 
highly nonlinear nature of the RS HJB PIDE and to the unboundedness of the instantaneous 
reward function g. 

This chapter is organized as follows. Section 2 introduces the general setting of the model 
and defines the class of random Poisson measures which will be used to model the jump com- 
ponent of the asset and factor dynamics. In Section 3 we formulate the control problem and 
apply a change of measure to obtain a simpler auxiliary criterion. Section 4 outlines the prop- 
erties of the value function. In Section 5 we show that the value function is a viscosity solution 
of the RS HJB PIDE before proving a comparison result in Section 6 which provides uniqueness. 



Our analytical setting is based on that of [19] . The notable difference is that we allow the factor 

processes to experience jumps. 

2.1 Overview 

The growth rates of the assets are assumed to depend on n valuation factors Xi{t), . . . ,Xn{t) 
which follow the dynamics given in equation @ below. The assets market comprises m risky 
securities Si, i = 1, . . . m. Let M := n + m. Let (0, {J^t} ,J^,^) be the underlying probability 
space. On this space is defined an M^-valued (J^t)-Brownian motion W{t) with components 
Wk{t), k = 1, . . . , M. Moreover, let (Z, Bz) be a Borel spac43- Let p be an (J^j)-adapted cj-finite 
Poisson point process on Z whose underlying point functions are maps from a countable set 
Dp C (0, oo) into Z. Define 



Consider Np{dt, dz), the Poisson random measure on (0, oo) x Z induced by p. Following Davis 
and Lleo |19j . we concentrate on stationary Poisson point processes of class (QL) with associated 

is a standard measurable (metric or topological) space and Bz is the Borel cr-field endowed to Z. 



2 Analytical Setting 
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Poisson random measure Np{dt,dx). The class (QL) is defined in |27j (Definition II. 3.1 p. 59) 

as 

Definition 1. An (J^()-adapted point process p on (0, J^, P) is said to be of class (QL) with 
respect to (J^t) if it is cr-finite and there exists iVp = (^Np{t, U)^ such that 

(i.) for U £ 3p, 1 1—)- Np{t, U) is a continuous (7"t)-adapted increasing process; 

(ii.) for each t and a. a. uj €z 0,, U >—?■ Np{t, U) is a cj-finite measure on (Z,^?(Z)); 

(iii.) for U e\,t^ Np{t, U) = Np{t, U) - Np{t, U) is an (J't)-martingale; 

The random measure |iVp(t, C/)| is cahed the compensator of the point process p. 

Since the Poisson point processes we consider are stationary, then their compensators are of 
the form Np{t,U) = u{U)t where v is the cr-finite characteristic measure of the Poisson point 
process p. For notational convenience, we define the Poisson random measure Np(dt,dz) as 

Np(dt,dz) 

_ ( Np{dt, dz) - Np{dt, dz) = Np{dt, dz) - u{dz)dt =: Np{dt, dz) if z G Zq 
\ Np{dt,dz) if z E Z\Zo 

where Zq C Bz such that zv(Z\Zo) < oo. 



2.2 Factor Dynamics 

We model the dynamics of the n factors with an affine jump diffusion process 

dX{t) = (6 + BX{t-))dt + AdWit) + [ ^{z)Np{dt, dz), X{0) = x (4) 

where X{t) is the M'^-valued factor process with components Xj{t) and b G M", B G M"^", 
A := [Kij] , i = 1, . . . , n, j = 1, . . . , and i{z) e M" with -oo < < ii{z) < ^f^"^' < oo for 
i = 1, . . . ,n. Moreover, the vector- valued function ^(z) satisfies: 

/ \az)\^i.idz) <oo 

(see for example Definition II. 4.1 in Ikeda and Watanabe |27| where Fp and Fp'°'^ are respec- 
tively given in equations 11(3.2) and 11(3.5)) 



2.3 Asset Market Dynamics 

Let So denote the wealth invested in the money market account with dynamics given by the 
equation: 

^^ = {ao + A'oX{t))dt, So{0) = so (5) 

where ao S M is a scalar constant, Aq € R" is a n-element column vector and where M' denotes 
the transposed matrix of M. Note that if we set = and cq = r, then equation ([5]) can be 
interpreted as the dynamics of a globally risk- free asset. Let Si{t) denote the price at time t of 
the ith security, with i = 1, . . . ,m. The dynamics of risky security i can be expressed as: 

dSi{t) 



{a + AX{t))idt + J2^ikdWk{t)+ [ ^iiz)Np{dt,dz), 
Sj(0) = Sj, i = 1, . . . ,m (6) 
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where a e R"", A £ R™^", E := [uij] , i = 1, . . . ,m, j = I, . . . , M and 7(z) G satisfies 
Assumption [2] 

Assumption 2. 7(2) € M"^ satisfies 

- 1 < < 7i(-z) < tI"""" < +00, i = 1, . . . , m 

and 

-l<7r"<0<7r"^<+oo, z = l,...,m 
for i = 1, . . . , m. Furthermore, define 

S := supp(z^) G jBz 

and 

S := supp(z^0 7-i) GBiR"") 

where supp(-) denotes the measure's support, then we assume that Ili^iiTi™"! 7i"'*^] the 
smahest closed hypercube containing S. 

In addition, the vector-valued function 7(2:) satisfies: 

/ \"f{z)fu{dz) < 00 

As note in [T9], Assumption [2] requires that each asset has, with positive probability, both 
upward and downward jump and as a result bounds the space of controls. 

Define the set J7 as 

J := |/i eM"^ : -1 -/iV < VV'GS} (7) 

For a given z, the equation h'^{z) = —1 describes a hyperplane in M*". Under Assumption [2] ^ 
is a convex subset of . 

2.4 Portfolio Dynamics 

We will assume that: 

Assumption 3. The matrix SS' is positive definite, 
and 

Assumption 4. The systematic (factor-driven) and idiosyncratic (asset-driven) jump risks are 
uncorrelated, i.e S Z and i = 1, . . . , m, ^i{z)£_'{z) = 0. 

The second assumption implies that there cannot be simultaneous jumps in the factor process 
and any asset price process. This assumption, which will prove sufficient to show the existence 
of a unique optimal investment policy, may appear somewhat restrictive as it does not enable us 
to model a jump correlation structure across factors and assets, although we can model a jump 
correlation structure within the factors and within the assets. 

Remark 5. Assumption (jH) is automatically satisfied when jumps are only allowed in the secu- 
rity prices and the state variable X{t) is modelled using a diffusion process (see |19j for a full 
treatment of this case). 
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Let Qt := a{{S{s), X{s)),0 < s < t) he the sigma-field generated by the security and factor 
processes up to time t. 

An investment strategy or control process is an R™-valued process with the interpretation 
that hi{t) is the fraction of current portfolio value invested in the zth asset, i = 1, . . . ,m. The 
fraction invested in the money market account is then /io(t) = 1 — YlT^i 



Definition 6. An M^-valued control process h{t) is in class Ti if the following conditions are 
satisfied: 

1. h{t) is progressively measurable with respect to {S([0, t]) (gi Gt}t>o cadlag; 

2. P (^f^ \h{s)f ds < +oo^ = I, VT>0; 

3. h'{t)-f{z) > -1, Vt > 0,z e Z, a.s. di^. 
Define the set JC as 

/C := {h{t) G n : h{t) ejyt a.s.} (8) 

Lemma 7. Under Assumption\^ a control process h{t) satisfying condition 3 in Definition\^is 
hounded. 

Proof. The proof of this result is immediate. □ 
Definition 8. A control process h{t) is in class AiT) if the following conditions are satisfied: 

1. h{t) enyte [0,T]; 

2. Ex^ = 1 where Xt is the Doleans exponential defined as 

Xt ■■= exp{-e [ h{s)'^dWs-\e'^ I h{s)'T:i:'h{s)ds 

■/ [ ln{l-G{z,h{s);e))Np{ds,dz) 
Jo Jz 



+ 



+ / / {ln{l - G{z,h{sy,e)) +G{z,h{s);e)}u{dz)ds 
Jo Jz 



(9) 



and 



G{z,h;e) = 1- {l + h'-f{z)) (10) 

Definition 9. We say that a control process h{t) is admissible if h{t) £ A.(T). 

The proportion invested in the money market account is /io(i) = 1 — hi{t). Taking this 
budget equation into consideration, the wealth V{t, x, h), or V{t), of the investor in response to 
an investment strategy h{t) € "H, follows the dynamics 

= (ao + A'oX{t)) dt + h'{t) {a - gqI + {A - IA'q) X{t)) dt 
+/i'(t)SdVFf + / h'{t)j{z)Np{dt,dz) 



where 1 € R™ denotes the m-element unit column vector and with V{0) = v. Defining a := 
a — a^l and A := A — l^g, we can express the portfolio dynamics as 

= (ao + A'^Xit)) dt + h'{t) [a + AX{t)) dt + h'{t)mWt + j h'{t)-f{z)Np{dt, dz) 

fill 



3 Problem Setup 



3.1 Optimization Criterion 



We will follow Bielecki and Pliska and Kuroda and Nagai ^31j and assume that the objective 
of the investor is to maximize the long-term risk adjusted growth of his/her portfolio of assets. 
In this context, the objective of the risk-sensitive management problem is to find h*(t) £ A{T) 
that maximizes the control criterion 



J{t,x,h;e) := -^InE 



-einV {t,x,h) 



(12) 



Hence, 



where 



By Ito, the log of the portfolio value in response to a strategy h is 

lnV{t) = lnv + {ao + A^Xis)) +h{sy (a + AX {s)^ds-^ h{sy^^'h{s)ds 

+ / /i(s)'SdVF(s) 
Jo 

+ [ [ {ln{l + h{sy-f{z)) -h{syj{z)}u{dz)ds 
Jo Jzo 

+ [ [ In {I + h{sy^{z)) Np{ds,dz) 
Jo Jz 



-ein V(i) 



h 



gix,h;e) 



+ 



f " exp <; w g{Xs, h{sy e)ds ) Xt 



+ 1) /I'SS'/i - ao - A'qX - h'{d + Ax) 



(l + /i'7(z)) '-1 +h'^iz)lz,iz)\u{dz) 



and the Doleans exponential Xt given by Q. 



(13) 



(14) 



(15) 



3.2 Change of Measure 

Let be the measure on {Q, T) defined as 

dFl 



dl 



Xt 



(16) 



For a change of measure to be possible, we must ensure that the following technical condition 
holds: 

G(z,/i(s);e) < 1 
for all s € [0,T] and z a.s. dv. This condition is satisfied iff 



h!{s)-i^z) > -1 



(17) 



a.s. dv, which was already one of the conditions required for h to be in class 7i (Condition 3 in 
Definitional). 
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i"^ is a probability measure for h G A{T). For h G A{T), 

Wl' = Wt + e [ E'h{s)ds 
Jo 



'0 

is a standard Brownian motion under the measure and we define the compensated Poisson 
measure as 



Jz 



N^{ds,dz) 



Jz 

t 



JZ 



Np{ds, dz) 
Np{ds, dz) 



Jz 

t 



{l-G{z,h{s);0)}i^{dz)ds 



JZ 



1 + 



h'-f{z)) ^^u{dz)ds 



As a result, X{s), < s < t satisfies the SDE: 



dX{s) = f {X{s~),h{s);e) ds + AdW^ + / ^{z)N^{ds,dz) 



where 



f{x,h;e):=b + Bx-eA^'h+ / C{z) {l + h'-f{z)) ^ - Izoiz) v{d 



(18) 



(19) 



We will now introduce the following two auxiliary criterion functions under the measure P^: 
• the auxiliary function directly associated with the risk-sensitive control problem: 



i{v,x-h-t,T;e) = -\\n^'i:', 



exp|(9^ g{Xs,h{s);e)ds-e\nv 



(20) 



where [■] denotes the expectation taken with respect to the measure P^ and with initial 
conditions (t, x). 

• the exponentially transformed criterion 



Iiv,x,h;t,T;9) := E,%^ 
which we will find convenient to use in our derivations. 



exp{e I g{Xs,h{s);e)ds-e\nv 

t 



(21) 



We have completed our reformulation of the problem under the measure P^. The state 
dynamics ()18p is a jump-diffusion process and our objective is to maximize the criterion (|20p or 
alternatively minimize ()2ip . 



3.3 The HJB Equation 

In this section we derive the risk-sensitive Hamilton-Jacobi-Bellman partial integro differential 
equation (RS HJB PIDE) associated with the optimal control problem. Since we do not an- 
ticipate that a classical solution generally exists, we will not attempt to derive a verification 
theorem. Instead, we will show that the value function $ is a solution of the RS HJB PIDE 
in the viscosity sense. In fact, we will show that the value function is the unique continuous 
viscosity solution of the RS HJB PIDE. This result will in turn justify the association of the RS 
HJB PIDE with the control problem and replace the verification theorem we would derive if a 
classical solution existed. 

Let $ be the value function for the auxiliary criterion function I{v, x; h; t, T) defined in ()20p . 
Then $ is defined as 

$(t,x)= sup I{v,x;h;t,T) (22) 
heA{T) 



We will show that $ satisfies the HJB PDE 

— {t,x) + snpL'mt,X{t)) = (23) 

where 

x) = f{x, h; e)'D^ + ^tr {KK'D'^^) - ^(D$)'AA'D$ 

+ y I _i _ _ - /i; e) (24) 

D- = and subject to terminal condition 

^>(r,3;)=ln7; (25) 



Similarly, let $ be the value function for the auxiliary criterion function I(v, x; h; t, T). Then 
^ is defined as 

$(t,x)= inf i(v,x;h;t,T) (26) 

heA{T) 



The corresponding HJB PDE is 



— {t, x) + -tr (^AA'D"$(t, x)j + F(a;, D<^) 
+ y |'l(t, a; + e(z)) - 'l(t, x) - C'{z)D^{t, x)} zy(dz) = (27) 
subject to terminal condition 

^{T,x) = v-^ (28) 

and where 

H{s,x,r,p) = mi_^{b + Bx-9AT.'h{s)yp + 9g{x,h;e)r^ (29) 

for r £ M, p £ M" and in particular, 

$(t,x) = exp{-e'^{t,x)} (30) 



The supremum in ()23p can be expressed as 
sup L^<I> 

= (6 + 5x)' i:>^> + ^tr (AA'L>2$) - ^(D$)'AA'i:)$ + ao + 

^Iz{~^ (e-e(*(M+?W)-*(M)) _ _^'(,)Z)ci>iz„(z)| K^^) 
+ sup I - - (6* + 1) /I'SS'/i - Oh'T.A'D^ + /i'(a + ix) 

{(1 - e^'(^)I)^) [(1 + /i'T^)"' - l] + eh'^{z)H,{z)} v{dz)^ 
Under Assumption [3] the term 



(31) 



\{e + l)h'T.T.'h-eh'^h'D^ + h'{a + Ax) - [ h'j{z)lzoizMdz) 
2 Jz 



is strictly concave in h. Under Assumption 21 the nonlinear jump-related term 



simplifies to 



iz) 



which is also concave in /i Vz G Z a.s. dv. Therefore, the supremum is reached for a unique 
optimal control h* , which is an interior point of the set J defined in equation ([T]), and the 
supremum, evaluated at h* , is finite. 



4 Properties of the Value Function 

4.1 "Zero Beta" Policies 

As in [19] ■ we will use "zero beta" (0/3) policies (initially introduced by Black [16])). 

Definition 10 (0/3-policy). By reference to the definition of the function g in equation ()15p . 
a 'zero beta' (0/3) control policy h{t) is an admissible control policy for which the function g is 
independent from the state variable x. 

In our problem, the set Z of 0/3-policies is the set of admissible policies h which satisfy the 
equation 

h'A = -Aq 

As m > n, there is potentially an infinite number of 0/3-policies as long as the following assump- 
tion is satisfied 

Assumption 11. The matrix A has rank n. 

Without loss of generality, we fix a 0/3 control h as a constant function of time so that 

g{x,h;e) = g 

where 5 is a constant. 

4.2 Convexity 

Proposition 12. The value function $(t, x) is convex in x. 

Proof. See the proof of Proposition 6.2 in [19j. □ 

Corollary 13. The exponentially transformed value function $ has the following property: 
V(a;i,X2) G M^^K G (0,1,), 

$(t, KXl + (1 - k)x2) > $"(*, Xl)$^~'"(t, X2) (32) 

Proof. The property follows immediately from the definition of $(t, x) = — ^ln$(t, x). □ 
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4.3 Boundedness 



Proposition 14. The exponentially transformed value function <I> is positive and bounded, i.e. 
there exists M > such that 



< x) <M V(t, x) G [0, T] X 



Proof. By definition, 



inf E 



h,e 

t,x 



exp 



g{Xs,h{s);e)ds - Olnv 



> 



Consider the zero-beta policy h. By tlie Dynamic Programming Principle 

^{t,x) < e^[-A'^5(^(^),/»;'?)<^^-in»;] ^ ge[s(r-t)-in,;] 

which concludes the proof. 
4.4 Growth 



□ 



Assumption 15. There exist 2n constant controls h^,k = 1, . . . ,2n such that the 2n functions 
13^ ■.[0,T]-^ M" defined by 



and 2n functions : [0, T] — )• M defined by 



where 
q{t) :-- 



b-OAT. 



'h+ [ 
Jz 



a{t) 



1 + h^ 7(z; 



q{s)ds 



(33) 



(34) 



i^idz) (t) 



-^tr (AA'/3'='(t)/3'=(t)) + j |e^'«(^) - 1 - e'(^)/3'^'(t)} z^(d. 

'l + h'^'jiz 

exist and for i = 1, . . . , n satisfy: 



+- 



+l-9{e + i)h'''^^'h'' -eao-Oa + e [ ll 

2 JZ I 



+ /i'=7(z)izo(-z) ^i^(fi^-z) 



m < 
/^r+^t) > 



(35) 



where /3](t) denotes the j-th component of the vector /3*(t). 

Remark 16. Key to this assumption is the condition ([35]) which imposes a specific constraint on 
one element of each of the 2n vectors [3^ it). To clarify the structure of this constraint, define 
-M^ as the square n x n matrix whose i-th column (with i = 1, . . . n) is the n-element column 

vector /3*(t). Then all the elements fnjjij = 1, . . . , m on the diagonal of are such that 

Similarly, define as the square n x n matrix whose i-th column (with i = 1 , . . . n) is the 
n-element column vector /3"+*(t). Then all the elements = 1,. . . ,m on the diagonal of 



are such that 



m- 



Note that there is no requirement for either or to have full rank. It would in fact 
be perfectly acceptable to have rank 1 as a result of column duplication. 
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Remark 17. For the function (3^ in equation ()33p to exists, B must be invertible. Moreover, the 
existence of 2n constant controls h^,k = 1, . . . ,2n such that ()33p satisfies (|35p is only guaranteed 
when J = M". However, since finding the controls is equivalent to solving a system of at most 
n inequalities with m variables and m > n, it is likely that one could find constant controls 
after some adjustments to the elements of the matrices Aq,A,B or to the maximum jump size 
allowed. 

Proposition 18. Suppose Assumptionll5\ holds and consider the 2n constant controls h^,k = 
1, . . . , 2n parameterizing the An functions 

a'' : [0,r] ^M, /c= l,...,2n 

[3^ : [0,r] ^M", A: = l,...,2n 

such that for i = 1, . . . ,n, 

m < 

/3^+\t) > 

where /3j(t) denotes the j-th component of the vector /3*(t). Then we have the following upper 
bounds: 

in each element Xj, i = 1, . . . , n of x. 

Proof. Setting Z = M" — {0} and recalling that the dynamics of the state variable X(t) under 
the P^-measure is given by 

dX{t) = f{X{t-),h{t);9)+AdWt^+ [ ^{z)N^{dt,dz) 

we note that the associated Levy measure D can be defined via the map: 

p = z/o^-i (36) 

We will now limit ourselves the class H'^ of constant controls. By the optimality principle, 
for an arbitrary admissible constant control policy h, we have 



^it,x) <iix;h;t,T) <Bt,, 



exp|0 J g{Xs,h)ds — 9lnv 



:=Wit,x) (37) 



In this setting, we note that the function g is an afRne function of the afRne process X{t). 
AfRne process theory See Appendix A in Duffie and Singleton j24j . Duffie, Pan and Singleton [23] 
or Duffie, Filipovic and Schachermayer [21] for more details on the properties of affine processes) 
leads us to expect that the expectation on the right-hand side of equation (|37p takes the form 

W{t, x) = exp {a{t) + /3{t)x} (38) 

where 

a : t e [0, T] ^ M 
/3 : t e [0, T] ^ M" 

are functions solving two ODEs. 
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Indeed, applying the Feynman-Kac formula, we find that the function W{t,x) satisfies the 
integro-differential PDE: 

^ + (^b + BXs- eK-L'h + ^ £,{z) [(1 + h'-i{z)) - izo (-s)] i^idz)^ DW{t, x) 



+ 



]-tx{KMD'^W(t,x))+ f {W{t,x + ^{z)) -W{t,x) - ^'{z)DW{t,x)} iy{dz) 
2 Jz 



+eg{x,h;e)Wit,x) 
= 

subject to terminal condition $(T, x) = v~^. 
Now, taking a candidate solution of the form 

W{t, x) = exp {a{t) + I3{t)x} 

we have 

dW 

~~dt 
DW 

D^W = P'{t)P{t)W{t,x 



a{t) + P{t)x] W{t,x) 
P'{t)W{t,x) 



Substituting into the PDE, we get 

(a{t)+${t)x^ W{t,x) 

+ (^ + Bx-ekT.'h + j ^{z)\{l + h'-i{z)y^ - lz^{z)\u{dz)^ P'(t)W{t,x) 



+-tr (AA'/3'(i)^(t)) W{t, 



+ J^{Wit,x + az))-W{t,x)-^'{z)p'{t)W{t,x)}u{dz) 



+e(-{e + l) h'T.T.'h - ao - A'qx - h'{a + Ax) 



+ 



= 



{i + h'j{z)y 



+ h'j{z)lzo{z)>iy{dz)]W{t,x 



Dividing by W{t,x) and rearranging, we get 

(/5(i) + B'p'{t) - eA'o - eJi'A^ X 
= - (^a{t) +(^b- 9KT,'h + j i{z) [(1 + h'-i{z)y^ - izo(^)] i^{dz) ) ^'{t) 



+1-9 {e + 1) /I'ss'/i -Oao-ed + e 

2 ./z 



{l + h'^{z)) ''-1 +h'^{z)lz,{z)\v{d 



Since the left-hand side is independent from the right-hand side, then both sides are orthog- 
onal. As a result we now only need to solve the two ODEs 



^{t) + B'l3'{t) - OA'q - Oh' A = 
12 



(39) 



and 



+-e {9 + 1) /I'SS'/i -9ao-9a + 9 [ 
2 Jz 



a{t)+ [b-9AJ:'h + J ^{z)[{l + h'j{z)y- lzo{z)\ Hdz)j 
+hr (AA'/3'(0/3(t)) + J [e^^^'^ - 1 - Hdz) 

(1 + h'jiz))-' - l] + /i'7(z)iz„(^)| Hdz) 
= (40) 
to obtain the value of W{t,x). The ODE (j39p for /3 is hnear and admits the solution 

/3{t) = 9B-^ (l - e^(^-*)) (Aq + /I'^i) (41) 



(42) 



As for the ODE ()40p for a, we only need to integrate to get 

cT 



a(t) 



}{s)ds 



where 



q{t) := [b-9AT.'h+ / ^{z) {l + h'^[z)) - Izoiz) i^{dz)) /3'(t) 



+ 



^tr (AA'/3'(t)/3(t)) + ^ {e''«(^) - 1 - e'(^)/3'(t)} 



+-9 {9 + 1) /I'SS'/i - 9ao -9a + 
2 ./z 



(l + /i'7(z)) '^-1 +/i'7(z)izo(^) '^(d^) 



Observe that W{t,x) is increasing in Xj, the i-th element of x, if /3i > 0, and conversely, 
W{t, x) is decreasing in Xi if /3i < 

Equations ()1T]1 and (P2|) are respectively equations (p3]) and ([34|) from Assumption [T5l By 
Assumption [151 there exists 2n constant controls h^, k = 1, . . . ,2n such that for i = 1, . . . , n, 

/3:(i) < 
/3^+\t) > 

where /3j(t) denotes the j-th component of the vector 13^ (t). We can now conclude that we have 
the following upper bounds 



for each element Xi,i = 1, . . . ,n of x. 



□ 



Remark 19. To obtain the upper bounds and the asymptotic behaviour, we do not need the 
2n constant controls to be pairwise different. In fact, we need at least 2 different controls and 
at most 2n different controls. Moreover, we could consider wider classes of controls extending 
beyond constant controls. This would require some modifications to the proof but would also 
alleviate the assumptions required for the result to hold. 



Remark 20. For a given constant control h, equation (j39p is a linear n-dimensional ODE. How- 
ever, if in the dynamics of the state variable X{t), A and H depended on X, the ODE would 
be nonlinear. Once ODE (|39p is solved, obtaining a{t) from equation ()40p is a simple matter of 
integration. 
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Remark 21. For a given constant control h, given x G and t G [0, T], the solution of ODE ()39p 
is the same whether the dynamics of S{t) and X{t) is the jump diffusion considered here or the 
corresponding pure diffusion model. The converse is, however, not true since in the pure diffusion 
setting h G M™, while in the jump diffusion case /i G JT" C M*^. 

5 Viscosity Solution Approach 

In recent years, viscosity solutions have gained a widespread acceptance as an effective tech- 
nique to obtain a weak sense solution for HJB PDEs when no classical (i.e C^'^) solution can 
be shown to exist, which is the case for many stochastic control problems. Viscosity solutions 
also have a very practical interest. Indeed, once a solution has been interpreted in the viscosity 
sense and the uniqueness of this solution has been proved via a comparison result, the funda- 
mental 'stability' result of Barles and Souganidis [8j opens the way to a numerical resolution of 
the problem through a wide range of schemes. Readers interested in an overview of viscosity 
solutions should refer to the classic article by Crandall, Ishii and Lions |17) . the book by Flem- 
ing and Soner [25] and 0ksendal and Sulem [30j, as well as the notes by Barles [5] and Touzi |34j . 

While the use of viscosity solutions to solve classical diffusion- type stochastic control prob- 
lems has been extensively studied and surveyed (see Fleming and Soner [26] and Touzi [M]), this 
introduction of a jump-related measure makes the jump-diffusion framework more complex. As 
a result, so far no general theory has been developed to solve jump-diffusion problems. Instead, 
the assumptions made to derive a comparison result are closely related to what the specific prob- 
lem allows. Broadly speaking, the literature can be split along two lines of analysis, depending 
on whether the measure associated with the jumps is assumed to be finite. 

In the case when the jump measure is finite, Alvarez and Tourin [1] consider a fairly general 
setting in which the jump term does not need to be linear in the function u which solves the 
integro-differential PDF. In this setting, Alvarez and Tourin develop a comparison theorem that 
they apply to a stochastic differential utility problem. Amadori [3j extends Alvarez and Tourin's 
analysis to price European options. Barles, Buckdahn and Pardoux [6] study the viscosity solu- 
tion of integro-differential equations associated with backward SDEs (BSDEs). 

The Levy measure is the most extensively studied measure with singularities. Pham [33] 
derives a comparison result for the variational inequality associated with an optimal stopping 
problem. Jakobsen and Karlsen [29] analyse in detail the impact of the Levy measure's singular- 
ity and propose a maximum principle. Amadori, Karlsen and La Chioma [3] focus on geometric 
Levy processes and the partial integro differential equations they generate before applying their 
results to BSDEs and to the pricing of European and American derivatives. A recent article by 
Barles and Imbert [7] takes a broader view of PDEs and their non-local operators. However, 
the authors assume that the nonlocal operator is broadly speaking linear in the solution which 
may prove overly restrictive in some cases, including our present problem. 

As far as our jump diffusion risk-sensitive control problem is concerned, we will promote a 
general treatment and avoid restricting the class of the compensator u. At some point, we will 
however need u to be finite. This assumption will only be made for a purely technical reason 
arising in the proof of the comparison result (in Section 6). Since the rest of the story is still 
valid if I' is not finite, and in accordance with our goal of keeping the discussion as broad as 
possible, we will write the rest of the article in the spirit of a general compensator i^. 



5.1 Definitions 

Before proceeding further, we will introduce the following definition: 
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Definition 22. The upper semicontinuous envelope u*{x) of a function n at x is defined as 

u*{x) = limsupn(y) 

and the lower semicontinuous envelope of u{x) is defined as 

u^{x) = liminfii(y) 

Note in particular the fundamental inequality between a function and its upper and lower 
semicontinuous envelopes: 

Ut: < u < u* 



The theory of viscosity solutions was initially developed for elliptical PDEs of the form 

H{x,u,Du,D^u) = 

and parabolic PDEs of the form 

— + H(x,u,Du,D^u) = 
ot 

for what Crandall, Ishii and Lions |17) term a "proper" functional H{x,r,p, A). 

Definition 23. A functional H{x,r,p,A) is said to be proper if it satisfies the following two 
properties: 

1. (degenerate) ellipticity: 

H{x,r,p,A) < H{x,r,p,B), B<A 

and 

2. monotonicity 

H{x,r,p, A) < H(x, s,p, A), r<s 
In our problem, the functional F defined as 

F{x,p,A) := - sup ^tr (AA'^) 

- -p'Kk'p 

-g{x,h)} (43) 

plays a similar role to the functional H in the general equation (|43p , and we note that it is indeed 
"proper". As a result, we can develop a viscosity approach to show that the value function ^ is 
the unique solution of the associated RS HJB FIDE. 

We now give two equivalent definitions of viscosity solutions adapted from Alvarez and 
Tourin [Ij: 

• a definition based on the notion of semijets; 

• a definition based on the notion of test function 
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Before introducing these two definitions, we need to define parabolic semijet of upper semicon- 
tinuous and lower semicontinuous functions and to add two additional conditions. 



Definition 24 (Parabolic Semijets). Let u G USC{[0,T] x M") and {t,x) G [0,r] x M". We 
define: 

• the Parabolic superjet Vu~^ as 

:= GMxM" x5„ : 

u{s,y) < u{s,x) +p{s -t) + {q,y - x) + - {A{y -x),y - x) 



+o(|s - t\ + \y - x\^) as (s,y) -> (t,a;)| 



— 2 + 

• the closure of the Parabolic superjet as 

T^t^ '■= \ {P,Q, A) = lim {pk,qk,Ak) with {pk,qk,Ak) £Vl'^ 
and lim {tk,Xk,u{tk, Xk)) = (t, x, u{t, x)) 

Let u G LSC{[0,T] x M") and it,x) G [0,r] x M". We define: 

• the Parabolic subjet Vu as Vu ■= —Vu^, and; 

— 2 — — 2 — — 2 + 

• the closure of the Parabolic subjet as = —V^ 

Condition 25 (Condition on an Upper Semicontinuous Function u). Let {t^x) G [0,r] x M"- 
and {p, q, A) G V'^'^u{t, x), there are (p G C(M"), tp>l and i? > such that for 

((s, y), z) G (^R(t, x) n ([0, T] X M")) x Z, 

Condition 26 (Condition on a Lower Semicontinuous Function u). Let (t,x) G [0, T] x and 
(p,g,yl) G 'P'^'~u{t,x), there are 99 G C(M"), > 1 and > such that for 

((s, y), z) G (^R(t, x) n ([0, T] X M")) x Z, 



The purpose of these conditions on u and v is to ensure that the jump term is semicontinuous 
at any given point {t,x) G [0,T] x M" (see Lemma 1 and Conditions (6) and (7) in [Ij). In our 
setting, we note that since the value function <I> and the function x 1— t- are locally bounded, 
these two conditions are satisfied. 

Remark 27. Note that the jump-related integral term 

.1 (^e-e{u{s,yH{z))-uis,y)) _ _ ^/(^)^| ^(^^^ 
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is well defined when {p,q,A) G Vu ■ First, by Taylor, 



z 



z 



_ 1 f^^-9(uis,y+az))-u(s,y)) _ -^^ _ ^'(^z)qj u{dz) 

9 2 

{u{s, y + ^{z)) - u(s, y)) - - (n(s, y + ^(z)) - u{s, y)) 



q2 

+ ^ {u{s, y + i{z)) - u{s, y)f + i'{z)q \ u{dz) 

By definition of the Parabolic superjet Vu~^ , for t = s, the pair {q,A) satisfies the inequality 

u{s, y + C{z)) - u{s, y) - i'{z)q < ^^'{z)AC{z) + o{\az)\^) 
Similarly, by definition of the Parabolic subjet Vu~, for t = s, the pair {q,A) satisfies the 



inequality 



u{s, y + e(z)) - u{s, y) - i'{z)q > ^C'W^^W + 



Thus, if M is a viscosity solution, we have 



uis, y + - u{s, y) - i'{z)q = \i'{z)Ai{z) + o{\i{z)\^) 



and the jump-related integral is equal to 



I {u{s, y + i{z)) - u{s, y)f + \i'{z)Ai{z) + \ v{dz) 



which is well-defined. 



Definition 28 (Viscosity Solution (Semijets)). A locally bounded function u £ USC{[0, T] xM") 
satisfying Condition [25] is a viscosity subsolution of ([23]) . if for all x G M", u{T,x) < go{x), and 
for all {t,x) G [0,r] x M", {p,q,A) G 'P'^'+u{t,x), we have 

-p + F{x,q,A)- J (^e-^Wt-^+5W)-«(t:^)) _ - u{dz) < 



A locally bounded function u G LSC{[0,T] x M") satisfying Condition [26] is a viscosity 
supersolution of ([23]), if for ah x G M", u{T, x) > go{x), and for ah {t, x) G [0, T] x M", {p, q, A) G 
V'^'~u{t, x), we have 



A locally bounded function $ whose upper semicontinuous and lowersemicontinuous en- 
velopes are a viscosity subsolution and a viscosity supersolution of (|23p is a viscosity solution 
of ([21. 

Definition 29 (Viscosity Solution (Test Functions)). A locally bounded function n G C/ SC{[0, T] x 
M") is a viscosity subsolution of ([23]) . if for all x G M", u{T,x) < goix), and for all (t,x) G 
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[0,r] X M", V G C2([0,T] X such that u{t,x) = 'tp{t,x), u < ip on [0,T] x M"\{(t,x)}, we 
have 

_ ^ + F{x,Di;,D^i;) - I (e-^('^(*'^+«(^))-'^(*'^)) - l) - ^'{z)D^l^^ u{dz) < 



A locally bounded function v G LSC{[0,T] x M") is a viscosity supersolution of (|23|) . if for 
ah X G M", v{T,x) > go{x), and for all {t,x) G [0,r] x M", V G C'2([0,r] x R") such that 
v{t, x) = ip{t, x), V > on [0, T] X M"\ {{t, x)}, we have 

_ ^ + F{x,D^,D^i;) - I [-] (e-^(^(*'^'+«(^«-'^(*'^)) - l) - '^('i^) > 



A locally bounded function $ whose upper semicontinuous and lower semicontinuous en- 
velopes are a viscosity subsolution and a viscosity supersolution of (j23p is a viscosity solution 
of (123]). 

We would have similar definition for the viscosity supersolution, subsolution and solution of 
equation (f27|) . Once again, the superjet and test function formulations are strictly equivalent 
(see Alvarez and Tourin [T] and Crandall, Ishii and Lions jljj). 

Remark 30. An alternative, more classical, but also more restrictive definition of viscosity solu- 
tion is as the continuous function which is both a supersolution and a subsolution of (j23p (see 
Definition 5.1 in Barles [5]). The line of reasoning we will follow will make full use of the latitude 
afforded by our definition and we will have to wait until the comparison result is established in 
Section [6] to prove the continuity of the viscosity solution. 

5.2 Characterization of the Value Function as a Viscosity Solution 

To show that the value function is a (discontinuous) viscosity solution of the associated RS HJB 
PIDE (|23p . we follow an argument by Touzi [M] which enables us to make a greater use of 
control theory in the derivation of the proof. 



Theorem 31. ^ is a (discontinuous) viscosity solution of the RS HJB PIDE ([23]) on [0, T] xM", 
subject to terminal condition (|25p . 

Proof. Outline - This proof can be decomposed in five steps. First, we define $ as a log 
transformation of In the next three steps, we prove that is a viscosity solution of the 
exponentially transformed RS HJB PIDE by showing that it is 1). a viscosity subsolution, 2). a 
viscosity supersolution and hence 3). a viscosity solution. Finally, applying a change of variable 
result, such as Proposition 2.2 in [34j, we conclude that $ is a viscosity solution of the RS HJB 
PIDE (|23]) 

Step 1: Exponential Transformation 

In order to prove that the value function <I> is a (discontinuous) viscosity solution of (|23p . we 
will start by proving that the exponentially transformed value function $ is a (discontinuous) 
viscosity solution of (f27]) . 

Step 2: Viscosity Subsolution 
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Let {to,xo) eQ := [0,t] x and u G C^'^(Q) satisfy 



and hence 



= ($* - u){to,xo) = max {^*{t,x) - u{t,x)) 



(44) 
(45) 



on Q 



Let (tk,Xk) be a sequence in Q such that 

hm {tk,Xk) = {to,xo) 

K— >00 

hm $(tfc,Xfc) = ^*(to,2;o) 
and define the sequence {C}k Ck ■= ^{tk,Xk)—u{tk,Xk)- Since u is of class C^'^, lim^-^oo Ck = 0. 

Fix h ^ J and consider a constant control h = h. Denote by the state process with 
initial data X^^ = x^ and, for k > 0, define the stopping time 

Tk := inf |s > tk : (s - tk,X^ - Xk) i [0, b^) x a^„| 

for a given constant a > and where is the unit ball in M" and 

5fc := V^(l- ^{o}(Cfc)) +A:"^-?{0}fe) 
From the definition of r^, we see that limfc-^oo''"fc = ^o- 

By the Dynamic Programming Principle, 



e^p{e I g{Xs,hs;B)ds}^Tk,X^J 



where 'Sjtf,^Xk ['] represents the expectation under the measure P given initial data (tfc,Xfc). 
By inequality (|i5|) . 



$(tfc,Xfc) < Ei,,^, 
and hence by definition of ^fc. 



expie / g{Xs, hs)ds } u{Tk,X!;J 



tk 



U{tk,Xk)+(,k < E; 



Tk 



exple g{Xs, hs)ds } u{Tk,X^) 



I.e. 



Tk 



exple g{XsX)ds \ u{Tk,X'^) 



u{tk,Xk) 



Define Z{tk) = 9 f^^ g{Xs, hs)ds, then 

d{e^') := 9g{Xs,h,)e^=ds 

Also, by Ito, 

■ l— + Cu\ds + DuK{s)dWs 

+ / {u{s,X{s-)+i{z)) -u{s,X{s~))}Np{ds,dz) 
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duo 



for s G [tfe,rfe] and where the generator £ of the state process X{t) is defined as 

/:u{t,x) := f(t,x,h;eyDu + ^tT{AA\t,X)D^u) 



(46) 



By the Ito product rule, and since dZs ■ Ug = 0, we get 

d {uge^") = Ugd (e^") + e^'dus 

and hence for t G [tfe,Tfe] 



k\„Zt 



u{tk,Xk)e^'^ +9 I u{s,X^)g{X^,hs)e^^ds 



+ 



+ 



j'^ (^^{s, X^) + Cu{s, X^)e^^^ ds + j'^ Du'K{s)dWs 

[t, X\s-) + i{z)) - u (t, X'^is-)) } Np{dt, dz) 



Noting that u{tk,Xk)e^*k = u{tk,Xk) and taking the expectation with respect to the initial 
data {tk,Xk), we get 



-•'tk 



du 
9^ 



{s,Xs) + Cu{s,Xs)+eu{s,Xs)g{Xs,hs) e^'ds 



In particular, for t = t^, 

6 < ^tk,Xk[u{Tk,XrJe^^k] -u{tk,Xk)e^*k 

~ ~^^tk,Xk 



-•'tk 



{—{s,Xs)+Cu{s,Xs) + eu{s,Xs)g{Xs,hs) ) e^^ds 



and thus 



f < Yi^tk,xk,H^k,Xr,)e^-'^]-uitk,Xk)e^^>^) 



1 



E 



fdu 
9^ 



tk 



(s, Xs) + Cu{s, Xs) + eu{s, Xs)g{Xs, hs) e^'ds 



As k ^ oo, tk ^ to, Tk ^ to, ^ Q and 



1 

du 



(^^{s,Xs) + Cu{s,Xs) + eu{s,Xs)g{Xs,hs)^ e^'ds 



—{s, Xs) + Cu{s, Xs) + eu{s, Xs)g{Xs, hs) 
a.s. by the Bounded Convergence Theorem, since the random variable 



L r (—( 



s,Xs) + jCu{s,Xs) + eu{s,Xs)giXs,hs) e^^ds 



is bounded for large enough k. 

Hence, we conclude that since hg is arbitrary, 

du 



— {s, Xs) + Cu{s, Xs) + eu{s, Xs)g{Xs, hs) > 
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- —is,Xs)-Cuis,Xs) - eu{s,Xs)g{X,,hs) < 

This argument proves that <I> is a (discontinuous) viscosity subsolution of the PDE (|27|) on 
[0,t) X M" subject to terminal condition <l(r, x) = e^"'^^'^). 

Step 3: Viscosity Supersolution 

This step in the proof is a shght adaptation of the proof for classical control problems in 
Touzi [33]. Let {to,xo) £ Q and u £ C^'^{Q) satisfy 

= - u){to,xo) < - u)it, x) for Q\{to,xo) (47) 

We intend to prove that at {tQ,xo) 

^(t,x) + mf |/:''n(t,2;) -eg{x,h)^ < 

by contradiction. Thus, assume that 

^{t,x) + mf^{c''u{t,x) - egix,h)] > (48) 

at {to,xo). 

Since C^u is continuous, there exists an open neighbourhood Ms of (to, 2:0) defined for J > 

as 

Ms ■■= {{t, x) : {t- to, x-xq) £ {-5, 5) X 6^n, and (gHD holds} (49) 

Note that by (071) and since 4> > ^* > li, 

min ( — n I > 

For p > 0, consider the set of p-optimal controls satisfying 

i{to,xo,hP) <^{to,xo) + p (50) 



Also, let e > 0, e < 7 be such that 

mm (^-u) > See-^^^"' > (51) 

where is defined as 

Ms := max {—g{x,h),0) 

{t,x)eAfg' ,hejp 

for 

AT/ := {{t, x):{t- to, x-xo)e {-6, 6) x {( + <5)^„} (52) 



and 



C := max \\(,{z) 



Note that C < 00 by boundedness of ^(z) and thus Ms < 00. 
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Now let {tk,Xk) be a sequence in Afs such that 

hm {tk,Xk) = {to,xo) 

K— >00 

and 

hm ^(tk,Xk) = ^:^(to,xo) 
Smce ($ — u)(tk,Xk) — >■ 0, we can assume that the sequence {tk,Xk) satisfies 

\{^-u){tk,Xk)\<e, foTk>l (53) 

for e defined by ([?T|) 

Consider the e-optimal control /i^ , denote by the controlled process defined by the control 
process h^. and introduce the stopping time 

Tfc := inf {s > Tfc : (s,X^(s)) ^ A^} 

Note that since we assumed that — oo < < ^ '?f^'*^ < for i = 1, . . . , n and since i' is 
assumed to be bounded then X{t) is also finite and in particular, 

(4> - u){Tk,Xl{rk)) > (3>* - u){Tk,Xl{Tk)) > 3ee-^^*^^ (54) 



Choose A/"/ so that (r, X^(t)) G A/"/. In particular, since X'^{t) is finite then A/"/ can be de- 
fined to be a strict subset of Q and we can effectively use the local boundedness of g to establish 
Ms. 



Let Z{tk) = 9 j^l^ g{XI, hl)ds, since $ > and by (l53]l and (fMD . 

$(Tfc,X^,(Tfc))e^(^'^) - l>(tfc,Xfc)e^(*'=) 

> ^x(T,,X^,(Tfc))e^(-^) - 4.(t,,Xfc)e^(*^) +3ee-'^^^-^^-e^(-'^-) 

> r d(u{s,Xl{s))e^A+2e 



I.e. 



Since we assumed that 



^{t,x)-C''u{t,x)-eg{x,h) <0 
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Htk,xk) < ^Tk,Xl{Tk))e^^^^'^ - r dU{s,Xl{s))e^^)-2^ 
Taking expectation with respect to the initial data (tfe,Xfc), 

m,Xk) < Et„,, |.(rfc,l^(rfe))e^(^'=)-^"'d(tx(s,X^(5))e^^) 

Note that by the ltd product rule, 

d[u{s,Xlis))e^' 
= Usd {e^") + e^" dus 
= ^{t,x) + C^u(t,x) + eg{x,h) 



2e 



then 

and therefore 



Tfe 



2e 



< -2e + E 



exp ^ / g{Xs,hl{s))ds\^{Tk,Xl{Tk)) 



< -2e + i{tk,Xk,hl) 

< $(tfc,Xfc)-e 

where the third inequahty fohows from the Dynamic Programming Principle and the last in- 
equality follows from the definition of e-optimal controls (see equation ([50]) ) . 



Hence, equation (jl8|) . 

— {t,x)+\ni^[c''u{t,x)-eg{x,h)^ >0 
is false and we have shown that 

— {t,x)+\ni^[C^u{t,x)-eg{x,h)^ <0 



This argument therefore proves that <I> is a (discontinuous) viscosity supersolution of the 
PDE (I27|) on [0,i) x W subject to terminal condition ^{T,x) = e9o(^;T). 

Step 4: Viscosity Solution 

Since ^ is both a (discontinuous) viscosity subsolution and a supersolution of (j27|) . it is a 
(discontinuous) viscosity. 

Step 5: Conclusion 

Since by assumption <I> is locally bounded, so is In addition, ^p{x) = e~^^ is of class Ci(]R). 
Also we note that ^ < 0. By the change of variable property (see for example Proposition 2.2 
in Touzi [34j), we see that 

1. since $ is a (discontinuous) viscosity subsolution of (|27p . ^ = ^p^^ o $ is a (discontinuous) 
viscosity supersolution of (p3]) : 

2. since <^ is a (discontinuous) viscosity supersolution of ()27p . ^ = o <3) is a (discontinuous) 
viscosity subsolution of (p3|) . 

and therefore <I> is a (discontinuous) viscosity solution of (|23p on [0, t) x M" subject to terminal 
condition «l(r,x) = e9o(^;^). 

□ 

We also note the following corollary: 
Corollary 32. (i). ^* is a upper semicontinuous viscosity subsolution, and; 

(a). is a lower semicontinuous viscosity supersolution 

of the RS HJB PIDE ([23]) on [0,T] x W, subject to terminal condition (|25]) . 

As a result of this corollary, we note that <!>*, and ^ are respectively a viscosity subsolu- 
tion, supersolution, and solution in the sense of Definitions [55] and [521 
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6 Comparison Result 



Once we have characterized the class of viscosity solutions associated with a given problem, the 
next task is to prove that the problem actually admits a unique viscosity solution by establishing 
a comparison theorem. Comparison theorems are the cornerstone of the application of viscosity 
theory. Their main use is to prove uniqueness, and in our case continuity, of the viscosity so- 
lution. Although a set of, by now fairly standard, techniques can be applied in the proof, the 
comparison theorem per se is generally customized to address both the specificities of the PDE 
and the requirements of the general problem. 

We face three main difficulties in establishing a comparison result for our risk-sensitive con- 
trol problem. The first obstacle is the behaviour of the value function $ at infinity. In the pure 
diffusion case or LEQR case solved by Kuroda and Nagai [31], the value function is quadratic in 
the state and is therefore not bounded for x G M". Consequently, there is no reason to expect 
the solution to the integro-differential RS HJB PIDE ([23]) to be bounded. The second hurdle 
is the presence of an extra non-linearity: the quadratic growth term {D^)' AA'D^. This ex- 
tra non-linearity could, in particular, increase the complexity of the derivation of a comparison 
result for an unbounded value function. Before dealing with the asymptotic growth condition 
we will therefore need to address this non-linear term. The traditional solution, an exponential 
change of variable such as the one proposed by Duffie and Lions [22], is equivalent to the log 
transformation we used to derive the RS HJB PIDE and again to prove that the value function 
is a viscosity solution of the RS HJB PIDE. However, the drawback of this method is that, 
by creating a new zeroth order term equal to the solution multiplied by the cost function g, it 
imposes a severe restriction on g for the PDE to satisfy the monotonicity property required to 
talk about viscosity solutions. The final difficulty lies in the presence of the jump term and of 
the compensator z^. If we assume that the measure is finite, this can be addressed following the 
general argument proposed by Alvarez and Tourin [[Ij and Amadori [2]. 

To address these difficulties, we will need to adopt a slightly different strategy from the 
classical argument used to proof comparison results as set out in Crandall, Ishii and Lions jl7] . 
In particular, we will exploit the properties of the exponentially transformed value function $ 
resulting from Assumption [15] and alternate between the log transformed RS HJB PIDE and 
the quadratic growth RS HJB PIDE ([23]) through the proof. 

Theorem 33 (Comparison Result on an Unbounded State Space). Let u = e~^^ G USC{[0, T] x 
M") be a bounded from above viscosity subsolution of ()23p and v = e^^" € LSC{[0,T] x M") be 
a bounded from below viscosity supersolution of ()23p . If the measure v is bounded and Assump- 
tion [75| holds then 

u<v on [0, T] X 

Proof. Outline - This proof can be decomposed in seven steps. In the first step, we perform 
the usual exponential transformation to rewrite the problem for the value function $ into a 
problem for the value function The rest of the proof is done by contradiction. In step 2, 
we state the assumption we are planning to disprove. The properties of the value function $ 
related to Assumption [15] are used in Step 3 to deduce that it is enough to prove the comparison 
result for $ on a bounded state space to reach our conclusion. We then double variables in step 
4 before finding moduli of continuity for the diffusion and the jump components respectively in 
steps 5 and 6. Finally, we reach a contradiction in step 7 and conclude the proof. 

Step 1: Exponential Transformation 



Let u G USC{[0,T] x M") be a viscosity subsolution of ([23D and u G LSC{[0,T] x M") be a 
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viscosity supersolution of ([23]) . Define: 



By the change of variable property (see for example Proposition 2.2 in Touzi [34j), u and v are 
respectively a viscosity subsolution and a viscosity supersolution of the RS HJB PIDE (|27p for 
the exponentially transformed value function <I>. 

Thus, to prove that 
it is sufficient to prove that 



u<v on [0, T] X M' 
u<v on [0, T] X M' 



Step 2: Setting the Problem 

As is usual in the derivation of comparison results, we argue by contradiction and assume that 

sup [u{t, x) — v{t, x)] > (55) 

(t,x)G[0,r]xR" 



Step 3: Taking the Behaviour of the Value Function into Consideration 

The assertion of this theorem is that the comparison result holds in the class of functions 
satisfying Assumption I15[ As a result Proposition [T2] holds and we can concentrate our analysis 
on subsolutions and supersolutions sharing the same growth properties as the exponentially 
transformed value function $. By Propositions 1181 and 1141 

< u{t, x) < e"' (*)^ V(t, x) G [0, T] x R'^ 

< v{t, x) < e"'(*)+^'' V(t, x) G [0, T] x R" 

and 

lim u{t,x)= lim ?}(t, x) = Vt € [0, T] (56) 

for k = 1, . . . ,2n where and /J'^ are the functions given in Assumption [T5j Since (|56p holds 
at an exponential rate, then by Assumption (|55|) there exists i? > 0, such that 

sup [u{t,x) — v{t,x)] = sup [u{t,x) — v{t,x)] 

{t,a:)g[0,T]xR" {t,x)e[0,T] xBij. 

Hence, it is enough to show a contradiction with respect to the hypothesis 

sup [u{t,x) - v{t,x)] > (57) 

established on the set Q := [0, T] x Br. Before proceeding to the next step, we will restate 
assumption ()57p now needs to be restated in terms of u and v as 

sup [u{t,x) - v{t,x)] > (58) 



Step 4: Doubling of Variables on the Set Q 
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Let r/ > be such that 



N := sup [u{t,x) - v{t,x) - ip{t)] > 



where f{t) := ?. 



We will now double variables, a technique commonly used in viscosity solutions literature 
(see e.g. Crandall, Ishii and Lions llZj). Consider a global maximum point (t^,Xe,ye) € (0,r] x 
Br X Br =: Qd of 

u{t, x) - v{t, y) - Lp{t) - e\x - 

and define 

:= sup x) — v{t, y) — ip{t) — e\x — yp] > 

Note that N^: > for e large enough. Moreover, N^: > N and as e ^ oo. 

It is well established (see Lemma 3.1 and Proposition 3.7 in |17] ) that along a subsequence 

lim {te,Xe,ye) = ii,x,x) 

for some {t, x) € [0, T] x which is a maximum point of 

u(t, x) — v{t, x) — ip{t) 
Via the same argument, we also have 



as well as 



and 



In addition, we note that 



lim ejxg — i/ep = 



lim u{t^,Xe) = u(t,x) 



lim v{t):,Xe) = v{t, x) 



hm N, = N 



Applying Theorem 8.3 in Crandall, Ishii and Lions [17J at (t^, x^, y^), we see that there exists 
tte, € M and A^, G Sn such that 

(a,,e(xe - ye), A) G 7^!'^ 
{b„e{x,-y,),B,) eVl'~ 

and 



- 3e 



" / ■ 




" A, 





< 3e 


" / -/ " 


< 




/ 







-Be . 




-/ / 
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Thus, we have for the subsolution u 



z 

< 

and for the supersolution v, 



l) +eC'(z)(x,-y,)|z.(dz) 



+ 

> 



z 



1 (^^-e{viu,y.+az))-v{u,y.)) _ 1^ + ee'(z)(x, - yA u{dz) 



Subtracting these two inequahties, 



< F{y^,e{xe - ye),B^) - F{xe,e{xe - ye),A) 

_ y 1^ (^^-0{u{Mi(z))-uit.,x.)) _-^}J^ eC'(z)(x, - ?/,)| lyidz) 
= F{y^,e{xe - ye),B^) - F(xe,e(xe - ye),A^) 



^-e{u{U,Xe+^{z))~u(U,x^))'^ ^f^^^- 



Step 5: Modulus of Continuity 

In this step, we focus on the (diffusion) operator F. 



F{ye,e{xe - ye),B^) - F(x,,e(x - y),A^) 

sup <^ ef{t„ y„ h)' (x, - y,) + -tr - -e^ (x, - y,)' AA' (x, - y,) - g{y„ h) 

- sup I e/(t„ x„ h)' (xe - 2/e) + Jtr (AA'^, + «„) 
-o^^^ (^^e - Ve)' AA' (xe - ye) - 9{xe,h) 



< hiv {AA'B, - AA%)| + sup {e\f{t,,y,,h) - f{t,,x,,h)\\{x, - y,)|} 

+ sup {|5r(Xe,/l) - g{ye,h)\} 

< hiv{AA'A,-AA'B,)\ + sup {e\f{t,,y,,h) - f{t,,x,,h)\\{x, - y,)\} 
+ snp{\g{xe,h) - g{y^,h)\} 



Note that the functional / defined in (fT9]l satisfies 

\f{t^,ye,h) - f{t^,Xe,h)\ < Cf \ye - x^ 
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for some constant Cf > 0. In addition 



tr {AA'A, - AA'B,) 
tr 





" AA' 


AA' ■ 




" A, 





( 


AA' 


AA' 







-Be _ 





" AA' 


AA' ■ 




" / -/ " 


( 


AA' 


AA' 




-/ / 



Finally, by definition of g, 



< 3e tr 
= 



\g{ye,h) - g{Xe,h)\ < Cg \ye - Xe 



for some constant Cg > 0. Combining these estimates, we get 

F{y^,e{xe - ye),B^) - F(xe,e(x, - ye),A, 
< uj{e\ye- x^f + \y^- Xe\) 



(60) 



for a function uj{C) = C(, with C = max [C/,Cg]. The function uj : [0, oo) — > [0, oo), which 
satisfies the condition a;(0^) = 0, is called a modulus of continuity. 



Step 6: The Jump Term 

We now consider the jump term 



v{dz) 



(61) 



Since for e > large enough, u{t, x) — v{t, y) > then 

u{t^,Xe + ^{z)) - u{t^,Xe) + v{t^,y^) - v{t^,ye + S,iz)) < -{u(te, x^) - v{t^,ye)) + N 

by definition of N. Moreover, since = sup(j 3, ,^)gQ^ [n(t, x) - v{t, y) - ip{t) - e\x - y\^] > 0, 
then < u{t^^x^) — v{t^,y^) and therefore 

u{t„x, +^{z)) - u{t„ Xe) + v{t„y,) - v{t„y, + ^{z)) < N - N, 

for z G Z. Thus, 

^-e(u(t,,x^+^{z))-u{t^,x^)+v{t^,y^)-v{t^,y,)) y ^-0{v{t^,y,+^{z))-v{U,y,)+N-N^) 

and equation ([6T]) can be bounded from above by: 



< 



^-e(v{t^,yc+^{z))-v{t,_,ye)) _ ^-9{u{te,x,_+^{z))-u{t,_,Xc)+v{tc,Xc)-v{tc,Xc))^ i/{dz) 
^-e{v,(U,y,+az))-vit„y,)) _ ^-e(vit„y,+az))-vit„y,)+N-Ns)^ ^^^^^ 
^-e{vit,,y,+az))~vit„y,)) _ e-^(^-^^)^| z.(dz) 
^-e{-^llnviU,y,+az))-hiv{U,y,)]) ( ^ _ ^-^(iV-TV,)^ | ^^^^^^ 



v{t„y, + C{z)) 

V{te,ye) 



l-e 



-e(N-Ne) 



v{dz) 



(62) 
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By Proposition [HI and since v is LSC, then 3A > : < A < v{t,x) < C^\/(t,x) G Q. As a 
result, 

for some constant K > In addition, since the measure v is assumed to be finite and the 
function i— ?> is continuous, we can estabhsh the fohowing upper bound for the right-hand 
side of dSSl): 

< LOniN-N^) sup iy{Z) (63) 

(t,{/)e[0,T]xR" 

for some modulus of continuity cj/j related to the function i-^ 1 — and parameterized by the 
radius R > oi the Ball introduced in Step 3. Note that this parametrization is implicitly 
due to the dependence of N and iV^ on R. The term sup(^ y)g[o,T]xR" ^(Z) is the upper bound 
for the measure u. 

Step 7: Conclusion 

We now substitute the upper bound obtained in inequalities (f60|) and (f63|) in ([59|) to obtain: 



<f'{t,) < u{e\y,-xf + \y,-x,\)+UR{N - N,) sup u{Z) (64) 

(t,x)g[0,T]xR" 



Taking the limit superior in inequality ()64p as e — ?> oo and recalling that 

(1) . the measure z/ is finite; 

(2) . S,i{z),i = 1, . . . , m is bounded Vz € Z a.s. dv 
we see that 

z^(Z) < oo 

Then 



\imujR(N - N^u(Z) = 



which leads to the contradiction 



^'(t) = I < 



We conclude from this that Assumption 1581 is false and therefore 



sup [v{t,x) — u{t,x)] > (65) 

{t,x)£Q 



Stated differently, we conclude that 

u<v on [0, T] X 



□ 
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6.1 Uniqueness 

Uniqueness is a direct consequence of Theorem 1331 Another important corollary is the fact that 
the (discontinuous) locally bounded viscosity solution $ is in fact continuous on [0,T] x R". 

Corollary 34 (Uniqueness and Continuity). The function ^{t,x) defined on [0, T] x M" is the 
unique continuous viscosity solution of the RS HJB PIDE (|23p subject to terminal condition (|25|) . 

Proof. Uniqueness is a standard by-product of Theorem[33j Continuity can be proved as follows. 
By definition of the upper and lower semicontinuous envelopes, recall that 

< ^> < 

By Corollary [32] and $* are respectively semicontinuous superolution and subsolution of 
the RS HJB PIDE 1^ subject to terminal condition (|^ 

We note that as a consequence of Theorem [33] is that 

> $* 

and hence 

is a continuous viscosity solution of the RS HJB PIDE (|23p subject to terminal condition (|25p . 

Hence, $ = <I>* = and it is the unique continuous viscosity solution of the RS HJB 
PIDE ([23]) subject to terminal condition ([25]). 

□ 

Now that we have proved uniqueness and continuity of the viscosity solution $ to the RS 
HJB PIDE ([23]) subject to terminal condition ([25]), we can deduce that the RS HJB PIDE ([27]) 
subject to terminal condition (|28p also has a unique continuous viscosity solution. We formalize 
the uniqueness and continuity of $ in the following corollary: 

Corollary 35 (Uniqueness and Continuity). The function <I>(t, x) defined on [0, T] x is the 
unique continuous viscosity solution of the RS HJB PIDE ()27p subject to terminal condition (|28p . 

7 Conclusion 

In this chapter, we considered a risk-sensitive asset management model with assets and factors 
modelled using affine jump-diffusion processes. This apparently simple setting conceals a num- 
ber of difficulties, such as the unboundedness of the instantaneous reward function g and the 
high nonlinearity of the HJB PIDE, which make the existence of classical C^'^ solution unlikely 
barring the introduction of significant assumptions. As a result, we considered a wider class of 
weak solutions, namely viscosity solutions. We proved that the value function of a class of risk 
sensitive control problems and established uniqueness by proving a non-standard comparison 
result. The viscosity approach has proved remarkably useful at solving difficult control prob- 
lems for which the classical approach may fail. However, it is limited by the fact that it only 
provides continuity of the value function and by its focus on the PDE in relative isolation from 
the actual optimization problem. The question is where to go from there? A possible avenue 
of research would be to look for a method to establish smootheness of the value function, for 
example through a connection between viscosity solutions and classical solutions. Achieving 
this objective may also require changes to the analytic setting in order to remove some of the 
difficulties inherent in manipulating unbounded functions. 
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